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ABSTRACT 

Voltage biased, frequency multiplexed TES bolometers have become a widespread tool in mm- wave astrophysics. 
However, parasitic impedance and dynamic range issues can limit stability, performance, and multiplexing factors. 
Here, we present novel methods of overcoming these challenges, achieved through digital feedback, implemented 
on a Field-Programmable Gate Array (FPGA). In the first method, known as Digital Active Nulling (DAN), 
the current sensor (e.g. SQUID) is nulled in a separate digital feedback loop for each bolometer frequency. 
This nulling removes the dynamic range limitation on the current sensor, increases its linearity, and reduces its 
effective input impedance. Additionally DAN removes constraints on wiring lengths and maximum multiplexing 
frequency. DAN has been fully implemented and tested. Integration for current experiments, including the South 
Pole Telescope, will be discussed. We also present a digital mechanism for strongly increasing stability in the 
presence of large series impedances, known as Digitally Enhanced Voltage Bias (DEVB). 

Keywords: Digital feedback, baseband feedback, TES bolometers, Field-Programmable Gate Arrays 

1. INTRODUCTION 

A large number of mm- wave astrophysics experiments^^ employ frequency multiplexed, voltage biased transition 
edge sensor (TES) bolometers to convert the incident radiation power into a current, which is measured by 
superconducting quantum interference devices (SQUIDs). A simplified schematic of this readout scheme is 
shown in Figure [T] 

Due to the steep total power-resistance relation of the TES, a voltage bias results in strong electrothermal 
feedback (ETF)Pln essence, the sum of optical and electrical power deposited on the bolometer is kept constant 
at the level dictated by the ETF loop gain. Series impedances to the TES spoil ETF and decrease stability, as 
we further explore in §3.2| The main sources of series impedance are the input impedance of the SQUID and 
the inductance of the cryogenic wiring. 

In addition, the SQUID has a highly non-linear response and has limited dynamic range. Providing negative 
feedback to the SQUID linearizes it and suppresses its effective input impedance. This negative feedback was 
previously provided broadband using a SQUID flux locked loop shunt feedbadCEl (hereafter referred to as shunt 
feedback). However, the shunt feedback has several issues. First, it needs to provide strong negative feedback 
over the entire bandwidth, while maintaining stability. Managing phase shifts and rolling off the open loop gain 
to maintain the stability of this loop is a challenged This requirement currently limits the usable bandwidth for 
frequency multiplexing to ~ 1.3 MHz and restricts wiring lengths. In addition, the presence of a superconducting 
leg in the LCR comb (i.e. a latched TES) can cause the shunt feedback to go unstable.^ 

Here we describe an alternative feedback scheme known as Digital Active Nulling (DAN), which does not 
suffer from these limitations, allowing a higher multiplexing factor, reduction of device parameter requirements, 
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Figure 1. Simplified schematic showing the voltage biased, frequency multiplexed bolometer readout system. The digitized 
carrier waveform provides strong voltage bias to each TES at its respective resonant frequency, as set by the series LC 
filter. The SQUID provides a current measurement, which is amplified, digitized and demodulated. In previous operation, 
the digitized miller signals provided static sinusoidal feedback, in addition to broadband feedback from the shunt feedback. 
The latter, shown dashed, is omitted in DAN operation, where the miller signal actively zeroes the SQUID current across 
the bolometer bandwidth. The thick lines highlight the main sources of series impedance to the bolometer. 

and improved stability. DAN is therefore a key technology for enabling multi-kilopixel bolometer arrays, for 
future ground-based, balloon-borne, and satellite applications. 

DAN is similar to baseband feedbackPEH (BBFB), except for the electrical injection of the nulling signal, 
whereas BBFB is applied directly to the SQUID through a separate feedback coil. If the coupling constant 
(defined as the ratio of mutual inductance to the geometric mean of self inductances) of the input coil to the 
feedback coil were unity, this would be equivalent to DAN as far as the feedback loop is concerned. However, 
BBFB applications have so far had a small mutual inductance of the feedback coil to the input coil, in which 
case the increased linearity and dynamic range are preserved, but the SQUID input impedance suppression is 
equal to the coupling constant. 

Furthermore, in ^3] we present a novel method of suppressing the effect of impedance in series with the TES 
known as Digitally Enhanced Voltage Bias (DEVB). Here, we once again take advantage of the small bandwidth 
of the TES and digitally correct for the voltage drop across the stray impedance; effectively providing strong 
voltage bias at each TES bias frequency separately. 

Series impedances spoil the voltage bias. In the extreme limit where the series impedance dominates the 
total impedance, the TES is effectively current biased and is unstable to perturbations. An increase in current 
causes an increase in electrical power and hence in resistance, leading to an instability. As shown in Appendix 
\K\ for an effective series complex impedance z s and TES resistance -RteSi the stability criterion is given by 

dlogVrES < \z s \ ^ 

dlogi?TES ~~ |^s+-RtEs| ' 

Once the logarithmic slope of the transition exceeds the fractional effective series impedance, the TES will latch 
into a superconducting state, such that it is no longer useful as a detector. In §3.2[ we show that DEVB allows 
for stable operation in this regime. 
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Figure 2. Schematic diagram showing the firmware portion of the DAN signal path. Data rates are shown in green. Both 
readout paths are shown. The residual readout path can be used to characterize the properties of the loop, though as 
long as the loop is operating correctly, the DAN readout alone respresents the bolometer current. 



2. DIGITAL ACTIVE NULLING 

We have implemented DAN on the McGill digital frequency multiplexing (DfMux) platform.^] Each DfMux 
circuit board provides 8 digital-to-analog converters, 4 analog-to-digital converters, and a Xilinx Virtex-4 Field- 
Programmable Gate Array (FPGA) for signal processing tasks. Each digital-to-analog converter synthesizes 
a "comb" waveform at 25 MSPS, consisting of 16 sinusoidal terms with independent amplitude, phase, and 
frequency controls. The DfMux also provides 64 complex demodulator channels, and streams data from each 
of them over Ethernet at a programmable sampling rate which, for this work, we set to 191 Hz. The software 
stack consists of an embedded Linux system, running C code and controlled remotely with a user-friendly python 

DAN is implemented as a separate discrete time integral control loop at each bias frequency. For one such 
bias frequency, the feedback loop is shown in Figure [2j We choose integral control for its desirable property 
of having an effective loop gain proportional to 1/(5/, where Sf denotes the frequency separation from the bias 
frequency. This implies that it has infinite gain at Sf — 0, where the bulk of the signal (the bias voltage) resides. 
For stability, we require that the open loop gain falls below unity at 180° of phase shift, such that no poles of 
the closed loop transfer function have a positive real component. The phase shift is dominated by digital time 
delays, which, in our implementation, range from 5.5 /is - 10.6 fis depending on the bias channel. Thus, we can 
easily achieve stability by letting the open loop gain fall below unity at ~ 10 kHz, which implies that at the 
highest bandwidth that may be required for science, say 100 Hz, we still have a loop gain of 100 and hence 99% 
effective nulling. 

2.1 Experimental Results 

We initially characterized our implementation of DAN by bypassing the cryogenic components of the readout 
system and using a warm resistive network that has a similar transfer function to that of the cryogenic compo- 
nents. A measurement of the residual signal at the current sensor as a function of frequency (hereafter referred 
to as a network analysis) of DAN in this characterization setup reveals the closed loop transfer function shown 
in Figure [3] 
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Figure 3. Left: Network analysis showing the broadband closed loop response of DAN. At 16 arbitrarily chosen bias 
frequencies where DAN is operating, the transfer function of an injected signal dips to zero i.e. the signal is nulled. Right: 
Zoom on the network analysis near one bias frequency for a variety of gain settings. The effective bandwidth of the nulling 
loop is shown to increase with G, a commanded gain parameter proportional to the open loop gain of the DAN feedback 
loop. 



We then operated DAN with cryogenic bolometers, including the SQUID readout. We found that one 
advantage of the DAN mode of operation is that the time required to map out the TES transition is reduced to the 
trivial exercise of stepping in bias amplitude while recording the DAN output. This mapping is performed after 
each cryogenic cycle on the SPTpol experiment.^ Using DAN decreases this tuning time (which previously took 
tens of minutes) by an order of magnitude, while simultaneously measuring the voltage-resistance characteristic 
at a much higher resolution. This speedup, simplification and improvement of the bolometer tuning algorithm 
has been implemented and is currently one of the uses of DAN occurring daily on SPTpol, and is planned to be 
used on the EBEX balloon-borne experiment.™ 

In the cryogenic test setup, we verified that the loop is operating correctly by streaming the residual demodu- 
lator signal and seeing a highly suppressed white noise level that increases with frequency in the manner expected 
from the integral control loop. We also verified that the DAN signal is the same as the demodulator signal when 
the gain parameter is set to zero. In particular, we provide a small sinusoidal test signal, which applies time 
varying power to the TES. The resulting demodulated current measurement!^ is found to be identical between 
the cases where DAN is enabled and where it is disabled. In addition, we compare the white noise level between 
these two modes of operation and find that they agree. A histogram of the white noise levels for a cryogenic 
SQUID without TES bolometers attached is shown in Figure |4j 

2.2 Improvements due to DAN 

Operating the system in DAN mode has several advantages. First, the strongly suppressed SQUID input 
impedance results in stronger voltage bias, improving TES stability. 

DAN nulls the SQUID current near bias frequencies, improving SQUID linearity and stability due to the large 
reduction of the total signal amplitude at the SQUID input. Similarly, improved nulling reduces the dynamic 
range requirement on the ADC, avoiding overloads in the case of unstable bolometers or changes in the bolometer 
operating point. 

Removal of the shunt feedback improves system stability, allowing longer wire lengths between the 4K SQUID 
and room temperature electronics, and increasing the allowed bandwidth. Additionally, the shunt feedback in- 
stability to a single superconducting leg in the comb is eliminated. DAN does not interact across bias frequencies, 
so any instabilities are localized to one TES, instead of the whole comb. We have experimentally verified this 
by allowing various numbers of individual bolometers to go superconducting and found no measurable effect on 
the noise level of the remaining bolometers on the same comb. 

The increase in usable bandwidth, removal of sensitivity to single misbehaving bolometers, and reduced usage 
of SQUID dynamic range, all enable a higher multiplexing factor. This higher multiplexing factor and the longer 
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Figure 4. Histogram comparing the white noise level for 15 DAN channels operating on a SQUID, spaced between 200 kHz 
and 1.2 MHz, with DAN enabled and DAN disabled. The two distributions are statistically consistent, confirming our 
expectation that the DAN operation does not alter the noise properties of the readout system appreciably. 

allowed wiring length are key for reducing heat load on cold stages such that it falls within the capabilities of 
available space cryogenics, as well as enabling lower cost, higher pixel count ground based experiments. 

Both increasing wiring lengths and moving the bolometer channels to higher frequency causes an increase 
in the wiring inductance, which is a source of series impedance to the TES. We note that the presence of this 
inductance shifts the optimal bias frequency and that at this frequency, the resulting effective impedance in series 
with the TES is smaller than jw Lairing- The amplitude and phase of the effective complex series impedance are 
non-trivially related to the details of the LCR network. 

The series impedance due to wiring can be mitigated by moving this wiring inside the DAN feedback loop, 
suppressing it in a manner similar to the SQUID input impedance suppression. This involves increasing the 
number of wires from each SQUID to a cryogenic stage near the detectors from two to six. An alternate scheme 
to mitigate this wiring inductance is DEVB, which we discuss in f|3j 



While the SQUID input impedance is strongly suppressed by DAN, residual series impedance (e.g. from wiring 
inductance, capacitor ESR, or magnetic coupling to lossy material) can limit TES stability. This series complex 
impedance can be measured to good accuracy by tuning the bolometers into their transition, cooling them 
below their superconducting temperature, and stepping down the bias voltage until the TES latches into its 
superconducting state. At this point, the current is proportional to the applied voltage (as we have experimentally 
verified), giving a measure of the residual complex impedance. We will show how this impedance can be used to 
dynamically adjust the voltage bias such as to keep the TES voltage constant. 

TES bolometers respond to incident power with a time constant To = C/G where C is the heat capacity of 
the TES and G is the thermal conductance to the heat sink. Under voltage bias, this time constant is sped up 
by the loop gain of electrothermal feedback /P as 



DEVB actively controls the voltage across the TES on time scales faster than T e ff, providing a constant voltage 
bias across the TES while strongly suppressing the effective series complex impedance. 

We denote this series impedance as z s = R s +jX s . Given values of R s and X s , a constant voltage bias to the 
TES denoted Vtes can be provided by measuring the current through the bolometer and providing the voltage 
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Figure 5. Firmware diagram showing the DEVB signal path. Data rates are shown in green. The complex current 
measurement and output voltage amplitude are read out simultaneously. 



where / denotes the magnitude of the measured current and Vtes is the programmable voltage set-point. Since, 
as with DAN, there are digital delays between measuring / and providing VouTj we perform this computation 
at a fixed loop gain below unity in order to introduce an effective DEVB time constant and ensure stability of 
the feedback loop. The DEVB time constant is discussed further in |3.1| 

In addition to providing improved TES stability, the enhanced voltage bias across the LCR comb suppresses 
cross-talk due to that same series impedance. In the case of an inductive series impedance, this term is propor- 
tional to the ratio of the stray inductance to the inductance of the LC resonator E With DEVB, this effect is 
suppressed by the DEVB loop gain. 

3.1 Preliminary Implementation 

A hardware prototype of DEVB has been implemented based on a customized version of the DfMux firmware^ 
This firmware supports 64 bias and miller channels, as well as 68 single phase demodulator channels. The DEVB 
feedback loop is implemented on one of those bias channels, using two of the demodulator channels, aligned 
90°out-of-phase to provide I and Q current measurements at the bias frequency. The resulting complex signal is 
fed to a DEVB module, which calculates an output amplitude for the synthesizer block by explicitly computing 
Equation [3j using the quadrature sum of the I and Q measurements as the current amplitude. The synthesizer 
modulates the amplitude of the bias sinusoid using this output. The input to, and the output from, the DEVB 
module are further decimated to 191 Hz and streamed across the network for setup and analysis. 

Figure [5] shows the computation within the DEVB module. At the input to this module, fixed-point signals 
from the demodulator are converted into floating-point numbers for computation. The algorithm itself is sep- 
arated into two parts: calculation of a target amplitude, and a control element generating the feedback signal 
by stepping the current amplitude a fraction towards the target amplitude. We denote this fraction the DEVB 
gain. The control signal generated by these blocks is finally converted back into fixed-point numbers suitable for 
the synthesizer block. 

We run the control loop at 12 kHz and set the DEVB gain to 1/3, yielding an approximate DEVB time 
constant of 4 kHz; much faster than the ~ 100 Hz time constant of a typical TES. 
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Figure 6. TES transition shown as the provided voltage is slowly decreased with DEVB enabled and disabled. The TES 
latches into a superconducting state when it reaches 60% of its normal resistance with DEVB enabled, compared to 
latching at 90% of its normal resistance with DEVB disabled. The left panel shows the raw data, whereas the right panel 
shows the data once the stray series impedance is subtracted off. The arrows are used to show the latching event, which 
occurs faster than the data rate. 



3.2 Measurements 

In order to demonstrate the enhanced performance of TES detectors under DEVB, we test the stability of a 
detector in the presence of a large stray impedance as given in Equation [T] Specifically, we apply a large carrier 
voltage to a single detector in an LC resonance, designed to allow for frequency multiplexing. By setting the 
frequency of this voltage significantly away from the LC resonance we create an impedance in series with the 
bolometer. In this case we set the series impedance to 1.4 j fi, while the normal TES resistance is ~ 0.9 f2. We 
then reduce the applied voltage, dropping the TES into its transition. The TES is lowered through its transition 
until it latches and becomes superconducting. The resulting total impedance as a function of applied voltage 
and the resistance of the TES alone as a function of the voltage across the TES alone are shown in Figure [6j 
This test, performed with DEVB enabled and DEVB disabled, shows that the use of DEVB lowers the point at 
which the TES latches, demonstrating that DEVB allows the TES to drop much lower into its transition in the 
presence of large series impedance. 

3.3 Future Prospects for DEVB 

While we have demonstrated the concept of DEVB, extending this feedback mechanism to an array of bolometers 
requires several improvements to the firmware, as well as the readout scheme and feedback mechanism. 

Equation [3j as well as the computation of the magnitude of the complex demodulated current, contains a 
square root operation, which is challenging to implement efficiently on an FPGA. An alternate DEVB scheme 
would involve controlling both the I and Q components of the bias voltage Vout rather than controlling only its 
amplitude. The set point for Vout would then become 

V"o UT = Vxes +I I R S - I Q X S 
V Q _t/Q +lQR +T i X 

v O\JT — V TES + 1 s s ■ 

This complex DEVB scheme, when implemented in fixed-point logic, would be easily scalable to the full 64 
bias and complex demodulation channels using the current Xilinx Virtex-4 FPGA. 

The current implementation of DEVB does not support DAN. In future versions of DEVB, we would imple- 
ment both the DEVB and DAN modules simultaneously. We would choose to run DAN at a significantly faster 
effective time constant than the DEVB control loop, such that the SQUID current is nulled, regardless of the 
changes in Vout made by the DEVB module. This approach would allow us to employ the benefits of DAN, such 
as the higher usable bandwidth and reduced SQUID dynamic range requirement, while retaining the suppressed 
series impedance from wiring inductance without needing to increase the heat load from the 4K stage to the 
millikelvin stage by increasing the wire count. 
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Finally, it is possible to implement DEVB with frequency-dependent coefficients such that a frequency- 
dependent series impedance could be suppressed. This would allow narrower LC filters to be used, increasing 
the multiplexing factor for a given usable bandwidth and constraint on electrical cross-talk. 



4. CONCLUSION 

We have introduced DAN, a method of keeping the current through the sensing portion of the readout system 
very close to zero at frequencies near the TES bias frequencies. This digital feedback mechanism increases 
stability and linearity of the SQUIDs, while suppressing their effective input impedance. We have demonstrated 
a fully implemented version of DAN including stability results and nominal noise performance, and described its 
use on the SPTpol experiment. 

While this technique is similar to BBFB, the suppression of series impedance is unique to DAN. Removing 
the shunt feedback from the readout system increases the usable bandwidth for frequency multiplexing, enabling 
higher multiplexing factors and reducing wiring length constraints. Hence, DAN is a key technology for applying 
frequency multiplexed readout to satellite applications. 

However, the use of higher bias frequencies increases the series impedance due to wiring. This effect can be 
mitigated by increasing the number of wires per multiplexing module from the 4 K SQUID to (close to) the 
millikelvin stage from two to six. 

We also described DEVB, a digital feedback technique which suppresses series impedances to the TES by 
measuring them and providing the actively controlled bias voltage required to hold the TES at a constant 
potential. This technique has been demonstrated to increase TES stability, despite large (z s > -Rtes) complex 
series impedances, without increasing wire count. 



APPENDIX A. STABILITY CRITERION FOR SERIES IMPEDANCE TO THE TES 

The left hand panel of Figure [6] shows that the reason for the TES latching due to series impedance is well 
understood. The total impedance is seen to have a zero derivative with respect to the total voltage, so further 
lowering the total bias voltage causes the voltage across the TES to fall in a runaway manner. This instability 
occurs when 

d\z +R TES \ ^ Q _ 
dvouT 

We substitute Equation [3] which describes the model, and assume that -Rtes is a monotonically decreasing 
function of Vtes- We find that this condition is met when 

d\ogV TES = \z s \ 
d\ogR TES \z s + i? TES | 

and conclude that the TES will latch due to series impedance when the fractional stray impedance exceeds the 
logarithmic slope of the TES transition. 
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